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Abstract 

The quenching factor for proton recoils in a stilbene scintillator was measured 
with a 252 Cf neutron source and was found to be 0.1 - 0.17 in the recoil energy 
range between 300 keV and 3 MeV. It was confirmed that the light yield depends 
on the direction of the recoil proton. 

The directional anisotropy of the quenching factor could be used to detect the 
wind of the WIMPs caused by the motion of the earth around the galactic center. 
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1 Introduction 



There is substantial evidence that most of the matter in our Galaxy must be 
dark matter and exist in the form of Weakly Interacting Massive Particles 
(WIMPs) [1]. WIMPs can be directly detected through elastic scattering with 
nuclei in a detector. The typical recoil energy is of order of 10 keV. The detector 
is required to measure about 10 keV nuclear recoil energies and discriminate 
between nuclear recoils and 7-rays to reject background events. 

For reliable detection of WIMPs, conventional searches depend mainly on ob- 
serving the annual modulation of the counting rate arising from the orbital 
motion of the earth around the sun and the rotation of the solar system itself 
around the galactic center. However, this approach is subject to a small mod- 
ulation amplitude and a large systematic error induced by long-term changes 
of the condition of the detector and seasonal variations of environmental ra- 
diations. For an alternative method, it might be useful to observe the wind of 
the WIMPs caused by the motion of the earth around the galactic center [2,3] 
(Fig. 1). The effect of the wind is large (~ 244 km/s) and a small systematic 
error can be achieved by means of periodically rotating the detector. To detect 
the wind of the WIMPs, a directional detector is required. 

It is known that organic crystal scintillators such as anthracene and stilbene 
have directional anisotropics in their scintillation responses to heavy charged 
particles [4,5]. Protons and carbon ions are produced internally by WIMPs 
scattering and the scintillation responses depend on the direction of their 
motion. These anisotropics can be used to obtain directional information of 
the WIMPs [6]. 

In using scintillation detectors to observe WIMPs, one needs to measure the 
scintillation efficiency or the quenching factor of the nuclear recoil. For organic 
crystal scintillators, the anisotropies of the scintillation responses are expressed 
in the directional variations of the quenching factor. 

Measurements of the quenching factor for nuclear recoils have usually been 
performed with several MeV mono-energetic neutrons produced by nuclear 
reactions using accelerators [7,8]. The nuclear recoil was caused by elastic 
scattering with neutrons instead of WIMPs. We performed the measurement 
employing a 252 Cf neutron source which is easier to handle. 

In this paper, we report on the measurement of the quenching factor of a 
stilbene scintillator as a function of the recoil energy and its dependence on 
the direction of the proton recoil. 
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2 Experiments 



The size of the stilbene crystal was 2 cm x 2 cm x 2 cm. It was purchased from 
Amcrys-HQ. Stilbene is an aromatic hydrocarbon including two benzene nuclei 
and forms a molecular crystal at room temperature which has an anisotropic 
structure. 

The crystal was covered with GORE-TEX® to reflect the scintillation light 
and its cleavage plane was glued to a PMT (H1161-50, Hamamatsu Photonics) 
with optical grease. 

The experimental setup is shown in Fig. 2. We employed a 3.56 MBq 252 Cf as 
a neutron source located 60 cm away from the stilbene scintillator. 

The recoil energy is calculated by the kinematics as 



2 (1 + A - cos 2 9 - cos 9^A 2 - 1 + cos 2 9) 

Er = (TTaJz K 



where 9 is the scattering angle, E n is the incident energy of the neutron and 
A = (mass of nucleus) / (mass of neutron). 

The scattering angle was determined by coincidence measurements between 
the stilbene scintillator and a liquid scintillator located 65 cm away from the 
stilbene scintillator at 60 degrees. Since the neutrons from 252 Cf source are 
not monochromatic but have a continuous Maxwellian energy distribution, 
the incident energy of the neutron needs to be measured by the TOF method 
between the source and the stilbene scintillator. To determine the timing of the 
nuclear fission of 252 Cf, a prompt 7-ray produced by the fission was detected 
with the plastic scintillator located by the neutron source. The energy of the 
scattered neutron was measured by the TOF between the stilbene scintillator 
and the liquid scintillator. 

The quenching factor Q is calculated by 

Q = -^visible/ -Er, (2) 



where -E v isibie is the measured electron-equivalent energy of the stilbene scin- 
tillator as calibrated by 7-rays. 



1 URL: http: / /www. aincrys-h.com 
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The recoil angle is calculated by 

2 , 1 + A - cos 2 9 - cos 9 V 'A 2 - 1 + cos 2 # 
cos = — . (3) 

For proton recoil, cos0 = sin#. We measured the quenching factor for two 
directions of the proton recoil, perpendicular and parallel to the cleavage plane. 

A schematic block-diagram of electronics for this experiment is shown in Fig. 3. 
The pulse height signal of the stilbene scintillator was taken from the dyn- 
ode output of the PMT, while the timing signal was taken from the anode 
output. The timing signals of the three scintillators were used for start- and 
stop-signals for two TACs. The pulses of the amplifier and the TACs were 
simultaneously recorded with a 4ch digital oscilloscope Tektronix TDS3034B. 
The peak pulse heights of the signals were sensed and digitized by the oscil- 
loscope and were sent to the PC using an RS232C serial line. 



3 Analysis 



Fig. 4 shows a scatter plot of the observed -Edible signal versus TOF between 
the neutron source and the stilbene scintillator. The events around ns are 
due to the 7-rays from the neutron source. The proton recoil events were 
distributed between 20 and 100 ns. Carbon recoil events have lower E vis iu e , so 
that they are not distinguished from background events. 

The recoil energy can be calculated from eq. (1) or by subtracting the 
energy of the scattered neutron from the incident energy calculated from the 
corresponding TOF. We used the latter method to obtain the recoil energy. 
The energy of the scattered neutron should be proportional to the incident 
energy as seen in eq. (1). We rejected background events which did not satisfy 
this relation within appropriate errors. 

Fig. 5 shows a scatter plot of the observed -^visible versus energy of the scattered 
neutrons. In case of the proton recoil, neutrons deposit 75% of their energy in 
the stilbene, so that we disregarded the events below the dashed line shown 
in Fig. 5 which have large residual kinetic energies and small invisible for the 
proton recoil. These events were most probably induced by carbon recoils and 
accidental coincidence background. 

Fig. 6 shows a scatter plot of observed invisible versus expected recoil energy 
£r. Taking the average of -E v i S ibi e /-E'R of the events in each bin of E-&, the 
measured quenching factor is shown in Fig. 7 as a function of the recoil en- 
ergy. We measured the quenching factor for two recoil directions, parallel and 
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perpendicular to the cleavage plane. The former is shown by filled squares and 
the latter by open circles in Fig. 7. 



4 Results and Discussions 

The quenching factor for proton recoils in a stilbene scintillator was measured 
with a 252 Cf neutron source and was found to range from 0.10 to 0.17 with 
recoil energy from 300 keV to 3 MeV. The measurement was performed for two 
recoil directions and we confirmed that the quenching factor depends on the 
recoil direction. It was smaller for the recoils perpendicular to the cleavage 
plane than for the recoils parallel to it. A measurement of the scintillation 
responses for a-particles was reported in Ref. [5], where the responses for the 
two directions are qualitatively consistent with our result. 

For the carbon recoil, definite data were not obtained because of the low invisible 
and the large background at low energies which were induced by accidental 
coincidences from 7-rays. To measure the quenching factor in the low -EVisibie 
region, we would need to improve the signal-to-background ratio. 

For the dark matter search, one can obtain directional information through 
comparing the data collected with different orientations of the detector. We are 
planning to measure two energy spectra directing the crystal axes of maximum 
and minimum scintillation responses toward the direction of the motion of the 
earth around the galactic center. 
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R : density of the number of WIMPs 

V x = 
V y = 



d 2 R 




Fig. 1. Expected 2D velocity distribution of WIMPs relative to the earth. In this 
representation, the earth is moving toward the positive x direction with a velocity 
of 244 km/s. V x and V y are x and y components of WIMPs velocity V, respectively. 
RMS velocity of WIMPs is assumed to be 230 km/s. 




Fig. 2. Experimental setup. 
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Fig. 3. Block diagram of the electronics for this experiment. Signals from an amplifier 
and two TACs were sensed and digitized by the oscilloscope and were sent to the 
PC. 
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Fig. 4. Measured scatter plot of invisible vs - TOF. 
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energy of scattered neutron [keV] 

Fig. 5. Measured scatter plot of invisible vs. energy of scattered neutron. We disre- 
garded events below the dashed line which have large residual kinetic energies and 
small i? v i s ibie for proton recoils. 




o 1 1 1 1 1 1 1 1 1 

500 1000 1500 2000 2500 3000 3500 4000 

£r [keV] 



Fig. 6. Measured scatter plot of invisible vs. Er. 
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Fig. 7. Measured quenching factor for proton recoils on average as function of recoil 
energy. The quenching factor is smaller for the recoils perpendicular to the cleavage 
plane than for the recoils parallel to it. 
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